function  [error]=plotPEN(MATGEN,test,n)
        maxmin=0;
        %min=0;
      %for in=1:serienum
          for mes=1:12
             
          maxmin(mes,1)=max(MATGEN(mes,:,1));
          maxmin(mes,2)=min(MATGEN(mes,:,1));
          %mim(mes,)
          end
          
          %media del modelo
          
          for serie=1:size(MATGEN,2)
            for m=1:12  
                fxi(m,1)=mean(MATGEN(m,:,1));
            end
          end
        
        error=sqrt(sum((fxi-test).^2)/12);
          
          x=[1:12];
          %subplot(3,1,1)
           
           figure(n);
           h1=plot(x,MATGEN(:,1,1),'-k');
           hold on;
          for seri=2:size(MATGEN,2)
           plot(x,MATGEN(:,seri,1),'-k');
            hleg1=legend('escenarios');
            hold on;
          end
          
          h2=plot(x,test,'-bd','LineWidth',2);
          legend([h1 h2],{'escenarios','real'})
          
        %subplot(3,1,2)   
        figure(n+20)
        %plot(x,maxmin(:,1),'-ro',x,maxmin(:,2),'-.b',x,test,'-bd')
        
        %hleg1 = legend('max','min','real');
        %subplot(3,1,3)
        plot(x,fxi,'-ro',x,test,'-bd')
        hleg2= legend('media','real');
        

end
